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Abstract 

We calculate the spin stiffness of the S = | frustrated Heisenberg antiferro- 
magnet directly from a general formula which is evaluated in the Schwinger 
boson mean-field approximation. Both Neel and collinear ordering are consid- 
ered. For collinear ordering, we take the anisotropy of this phase into account, 
unlike previous approaches. For Neel ordering, a detailed study is made of the 
finite-size scaling behavior of the two terms that make up the spin stiffness. 
The exponents of the scaling with the system size of the two terms comprising 
the spin stiffness turn out to be identical to those of the unfrustrated case. 
PACS numbers : 75.10.Jm, 75.30.Kz, 75.40.Cx 

I. INTRODUCTION 

The recent interest in the frustrated Heisenberg antiferromagnets is motivated by high 
T c -superconductivity; the undoped compounds show long-range antiferromagnetic order, 
similar to the Heisenberg model. Upon doping superconductivity occurs. Adding frustration 
to the Heisenberg model can be thought of as to mimic the effect of hole doping. 

We consider the frustrated Heisenberg model on a square lattice with N = L 2 sites. It 
is described by the following Hamiltonian for quantum spins Sj on a lattice: 

H = J 1 ^S l -S J + J 2 5> t -S„ (1) 

nn nnn 

where nn denotes a pair (ij) of nearest-neighbor sites and nnn a pair of next-nearest-neighbor 
sites. The spin length is fixed; S = ^ Both J\ and J2 are taken to be non-negative. If 



J2/J1 is small, the antiferromagnetic long-range order is recovered (Neel-like). For J2/J1 
large, the system decomposes in two Neel ordered sublattices which, however, have the 
same quantization axis. Alternating strips of up and down-spins will occur; the so-called 
collinear ordering. Clearly these couplings frustrate each other. If the spins were classical 
objects a large number of phases, among which the Neel and collinear phases, would become 
degenerate for J2/J1 = 0.5. For the quantum case, a quantum phase transition to a spin- 
liquid phase might occur.Ehi 

In this paper, our intention is to employ the spin stiffness p s to measure magnetic order 
in the system. We calculate the spin stiffness in the framework of the Schwinger boson 
mean-field approximation (SBMFA) using a general formula for p s . Previous evaluations of 
the spin stiffness were indirecti§§. However, apart from minor adjustments, we confirm 
their results. Furthermore the scaling behavior of p s in this approximation is derived. This 
is useful to sensibly extrapolate results of more exact approaches, like quantum Monte Carlo 
and exact diagonalization. Also, the typical system sizes for which scaling is valid can be 
estimated in this way.i~@ 

II. SCHWINGER BOSON MEAN-FIELD APPROXIMATION 

The SBMFA improves upon standard mean-field theory by incorporating correlations 
between neighboring spins00. With this approximation we derive the energies and wave- 
functions of all states of the frustrated Heisenberg model for both Neel and collinear order. 
Our notation below will generally follow that of Mila et al.H 

The Schwinger boson transformation is a representation of the separate spin operators 
by pairs of boson operators; S + = a^b, S~ = atf , S z = \ a^a — tfb , supplemented by the 
local constraint a^a + b^b = 2S(= 1). 

To transform the Hamiltonian in a convenient form an appropriate rotation in spin space 
is applied. Define Dij = a^a] + bM and = a^b^ + biaj. The Hamiltonian becomes 
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h = -i E ^(4^ - 1) + 1 E - 1), (2) 

AFM FM 

where we have inserted for the pairs {if) of spins parallel in the S^-direction and 
for the antiparallel pairs. For the two orderings considered this is depicted in Figure |]. The 
parameter equals J\ for nearest-neighbors and J 2 for next-nearest-neighbors. The mean- 
field decoupling is made using the fields Ky = \(Dij) and 7^ = The local constraint 
a^a + b^b = 1 is replaced by a global one and enforced by means of a Lagrange multiplier 
A. After a Fourier transform from (ai,bi) to (a p ,b p ) and a Bogoliubov transformation from 
(a P ,6p) to (a p ,Pp); 

a P = a P cosh 9 p + f3 p sinh p , 
bL p = a p sinh 9 P + /3_ cosh p , 

the Hamiltonian becomes 

H MF = E C + J2 w P («t«p + /3p/3 P ), (4) 
p 

where we have introduced the quasiparticle energies 



^p = V^p + A ) 2 - A p- ( 5 ) 

Below the quantities h p , A p and E c will be defined for the Neel and collinear orderings 
separately. The fields k^j, 7^ and the Lagrange multiplier A are obtained through consistency 
equations. 

We consider two types of order: 

Neel order: the specific form of h p , A p , E c and the three consistency equations are: 

h p = 4J 2 K cosp x cospy, (6) 
A p = 2J 1 7(cosp a; + cospy), (7) 



E c = 2N 



J 1 (i 4+ 2 1 2 )-J 2 (I + 2k 2 )-\ , (8) 
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K = ^J2~^~ C0S P* C0S PV> ( 9 ) 

7 = ^H^f( cos ^ + cosp y ), (10) 

For future applications it is also useful to define the "condensate" mjl by 

h + A _ 1 „ Zip + A 



1__L y (12) 



This is the combination of the -equivalent- p = (0, 0) and p = (tt, it) terms in ([5J) and 
(PI), which both diverge for N — > oo. In the same limit, N — > oo, m s also equals the 
corresponding terms in fllPp. 

Collinear order: for this phase we introduce quantities with a bar where confusion might 
arise; 

hp = 2J\R cosp x , 

A p = 2 J171 cosp y + 4J 2 7 2 cosp x cosp y , (13) 

E c = 2N (Jx( 7l 2 - R 2 - i) + J 2 (i + 2 72 2 ) - A) , (14) 

R = — p — cosp x , (15) 
A^ 2cj p F ' v ; 

71 = — — — cosp,,, (16) 

72 = ^Xl^-COSPajCOSPj,, (17) 

iV p ZuJ p 

The condensate m s and cjp are defined in a similar manner as before in (|T2"D and (||). 

The discussion above has given us the ground states |0) with energy E = E c + I] p o; p 
for both orderings. These ground states are characterized by the absence of quasi-particles; 
a\0) = P\0) = 0. Excited states are given by: 

\a) = a\-...-ai-pl-...-Pl\0), (19) 



H M f\o) = E a \a) 

= [E Q + ujx + . . . + uj n + ujx + . . . + u m ]\a) . 



(20) 



III. THE SPIN STIFFNESS IN THE SBMFA 

The spin stiffness p s is non-zero if there exists magnetic order in the system and is, 
at T = 0, associated with an increase in energy upon twisting the order parameter of the 
system (AE = |iVp s |q| 2 with q the wave vector of the twist. 

00. In line with Einarsson 
and SchulJ we introduce this twist through a replacement of Dij and Bij in (§) by 

Ai(q) = a ia ]e^ rs + b^e'^ (21) 
B y (q) = o^e*™ + ka 3 e-^ rs , (22) 

with r$ = Tj — rj. The resulting Hamiltonian H(q) is now evaluated within the SBMFA. 
Defining /%-(q) = |(Dy(q)) and 7#(q) = |(-By(q)), the mean-field Hamiltonian becomes 

H M f(vl) = - E ^7y(q)(4i(q) + 5 y (q) - 2 7ii (q)] 



AFM 



+ E ^v«w(q)Py(q) + Ai(q) - 2«y(q)] 



FM 



+A E( a l a « + 4^ — 1) + constants. (23) 

i 

Since within the SBMFA we know what the excited states are [see (|19| ) and (|20|) 1, the spin 
stiffness can be directly evaluated from second order perturbation theory 

» = ->M0> + £ £ «^ + ,, (24) 
with the quantities t and j defined by 



t = -^HmM) 



, 3 = ■^H MF (q) 
q =o d 1 



(25) 



q=0 



In (24), we have also defined the abbreviations T and J for the two terms in p s . 
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IV. RESULTS FOR THE SPIN STIFFNESS 



From this point on, we set q = g(cos</>, sin0). Of the two terms for p s in fl2"j|), T is 
evaluated more easily. We obtain after straightforward manipulations with Brillouin zone 
summations : 

T = J l7 2 - 2J 2 k 2 , (26) 
f = 2J 27 2 + Ji( 7l 2 sin 2 <p-n 2 cos 2 <p), (27) 

These simple equations hold for all system sizes N. 

The quantity J requires more effort; the operator j has to be expressed in the operators 
ct p and j3 p defined in fl3|). For the wavefunctions \a) and energies E a of the excitations we 
use ( |19D and (p0|) . The resulting values for J are written as summations over the Brillouin 
zone: 

J = -Jj E S -^T 1 (JMh P + A) - 2 J 2 kA p cos Px ) 2 , 

1 1 _ - - 

J = -j7 22 — (cos0sinp x [JiKA p - 2J 2 7 2 (/i p + A) cosp^] 

p w p 

- sin^sinpj,[Ji7i(^ p + A) + 2J 2 7 2 (/i p + A) cospj) 2 . 

Only for the infinitely large lattice these equations can be simplified by replacing summations 
by integrals and partially integrating. The expression for the spin stiffness p s = T + J then 
simplifies considerably and becomes: 

p s = m s (Ji7 - 2J 2 k), (28) 
p s = m s [(2J 2 7 2 - JiR) cos 2 (j) 

+ (2J 2 72 + Ji7i)sin 2 0]. (29) 

Ivanov and Ivanovil apply a different method to derive p s . They use a modified spin wave 
theory which leads to the same consistency equations as the Schwinger boson approach. 
p s is then obtained by calculating the correlation length £ associated with the spin-spin 
correlation function (Sj • Sj) and comparing this £ to the expression for £ obtained for the 
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non- linear sigma model to two- loop order by Chakravarty et al .1131 (where £ ~ exp(2np s /9) 
with 9 the temperature). For Neel ordering their expression is identical to our result (pBp . 
It is gratifying to see that the non-linear sigma model also is the effective field theory for 
the low-energy physics of the frustrated Heisenberg antiferromagnet. 

For the collinear ordering they obtain the geometrical average of our cos 2 (j) and sin 2 en- 
terals, whereas we take explicitly the anisotropy of this phase into account, Still both ex- 
pressions for p s vanish at the same value of Jij ' J\. 

Table | and Figure contain our numerical results. 

V. SCALING OF THE SPIN STIFFNESS 

It is necessary to know the size dependence of observables to obtain a good approx- 
imation for their limit values. Neuberger and Zimanlil derived the scaling behavior for 
an unfrustrated Heisenberg antiferromagnet explicitly. Here we extend this to the case of 
frustration. 

Recently some discussion has arisen about where the scaling behavior of p s sets ini. Our 
formulas in the last section lend themselves well to investigate this. 

Here we only treat the Neel ordering. We want to know the scaling behavior of the 
condensate m s and the two terms J and T that make up p s (p s = T + J). The latter two 
will turn out to have different scaling behavior. 

As can be seen from (^) only the scaling behavior of kn and 7at is required for T. These 
two are part of the set (kat, 7jv, A at) of mutually dependent quantities. We will now argue 
what is the exponent of their scaling behavior and therefore of T, without trying to obtain 
the precise prefactor (which would be quite tedious). 

Name the equations ([)])-(|ll]) /, 77, and III, respectively. They contain poles at p = (0, 0) 
and p = (7r,7r). With help of ([Sp, (]B|) and (0) we rearrange them as / — III, 77 — III 
and 4J 2 knI — ^JiInII + \nIII. We neglect the p = (0,0) and p = (tv,tt) terms in the 
summations. It is easy to show that this will give rise to errors of the order 0(N' 2 ). Next 
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we expand these equations to first order round their infinite size values (k, 7, A). Define the 



size dependences 5k n = 


= K N - K, 5-f N = 7iv - 


- 7 and SXn 


— Xn 




K-l 




Skn 




7-1 


= A + B- 


8 In 




A - 4J l7 2 + 4J 2 k 2 




5X N 



(30) 



where A and B contain summations over the Brillouin zone dependent on the infinite size 
parameters k, 7, and A. The remaining size dependence of B can be neglected as it leads to 
higher order terms. On the other hand the summations in A will be replaced by integrations 
plus size dependent corrections. Using Neuberger and Ziman0 we obtain 



B 



£>lN 

5X N 



1 c. 



(31) 



The parameters k^, 7^ and Xn thus scale with N~ 3 ^ 2 . A direct consequence of this is that 
T N — T ~ iV~ 3 / 2 . If the size dependence of the parameters is neglected [(k-n,7n, ^n) 
(k, 7, A) inside the summations (|9|)- ([TT|) 1, we also find the prefactor: 



T N -T 



0.7186 1 X + AJ 2 k 



[A - 8 J 2 k] 



(32) 



iV 3/2 V2 V A - 4 J 2 k 

Upto J 2 /Ji ~ 0.5 this is in good agreement with the numerical solution of ©-(0) for 
various sizes N. 

Next we consider the condensate m s n and Jn, for which we can even derive the prefac- 



tors. Replacement of (kn, 7Ar, A^v) by («, 7, A) in the summation of ( [T2] ) leads to errors of the 
order 0(iV -3 / 2 ). We will neglect these. Neuberger and Ziman^l state a lemma applicable to 
this summation. This leads to: 

0.6208 r; IX + AJ 2 k 



m S N - m s 



For Jn we find in a similar fashion: 



N 



X - 4J 2 k' 



(33) 



Jn — J 



0.6208 V A2 - (4</ 2 «0 2 
4^2 ' 



(34) 



These formulas are in excellent agreement with the numerical solution to (|9|)-(|lT|). For Jn 
this is depicted in figure 0. 

In conclusion, we have obtained the scaling behavior of T and J from an analysis of 
the formulas in the SBMFA for the frustrated Heisenberg antiferromagnet. The qualitative 
scaling behavior (i.e. the exponents) is the same as for the unfrustrated case, which was 
discussed by Neuberger and Zimanlll. The scaling behavior that was utilized by Einarsson 
and Schulzi does agree with our findings but we confirm the message of Feiguin et al.@ that 
the clusters they used are too small for the scaling behavior of J to have set in. Numerically 
we see that the scaling behavior starts around size N = 100 (See Figure |2|) whereas the 
largest cluster they used is N = 36. 

We would like to thank J.M.J, van Leeuwen and W. van Saarloos for fruitful discussions. 
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TABLES 



J2/J1 


Order 


T 


-J 


Ps 


0.0 


N 


0.3352 




0.1596 




0.1757 




0.1 


N 


0.2961 




0.1596 




0.1365 




0.2 


N 


0.2597 




0.1597 




0.1000 




0.3 


N 


0.2271 




0.1600 




0.0672 




0.4 


N 


0.1995 




0.1604 




0.0391 




0.5 


N 


0.1783 




0.1612 




0.0171 




0.6 


N 


0.1639 




0.1622 




0.0017 




0.6 


C 


0.6231 


0.0499 


0.3841 


0.0284 


0.2390 


0.0214 


0.7 


C 


0.7211 


0.1563 


0.3327 


0.0797 


0.3884 


0.0766 


0.8 


c 


0.7794 


0.2566 


0.3232 


0.1264 


0.4563 


0.1302 


0.9 


c 


0.8327 


0.3492 


0.3292 


0.1689 


0.5034 


0.1803 


1.0 


c 


0.8866 


0.4361 


0.3434 


0.2086 


0.5432 


0.2275 


1.1 


c 


0.9420 


0.5190 


0.3623 


0.2465 


0.5797 


0.2725 



TABLE I. The limit (N = 00) values for T, —J and p s as function of the ratio J2/J1 where 
J\ = 1. The two orderings considered are Neel (N) and collinear (C) order. For the collinear 
ordering there are two directions: along the antiparallel spins, 4> = 7r/2 (listed first) and along the 
parallel spins, <j) = (listed secondly). 
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FIGURES 



0.6 




0.2 0.4 0.6 0.8 1 

J 2 /Ji 



FIG. 1. The spin stiffness p s in units of J\ (solid lines). For the collinear ordering the spin 
stiffness in the direction of the parallel spins, = 0, (lower solid curve) and in the direction of the 
antiparallel spins (j) = it/ 2 (upper solid curve) are drawn. The dotted line is the result found by 
Ivanov and Ivanovi. Inset : The mean fields for the Neel (a) and the collinear order (b). Mind 




0.2 0.4 0.6 

Jz/Jl 

FIG. 2. The numerical scaling behavior of Jm (solid lines) for sizes L = 2,4, 6, 8, 10, 20, 40, 100 
(numbering is top-down) compared with the theoretical formula ( |34| ) (dotted line). 
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